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This  report  describes  a  summary  of  the  use  of  probabilistic  event  combination  techniques 
in  the  design  and  evaluation  of  navigation  structures  and  presents  recommendations  for 
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1.  INTRODUCTION 


1.1  Review  of  Previous  Work 

Navigation  structures  must  be  designed  to  withstand  the  effects  of  loads 
and  their  combinations  arising  from  operating  conditions,  environmental  events 
and  effects,  and  accidents  such  as  barge  impact.  The  loading  requirements  in 
design  documents  used  by  the  U.S.  Army  Corps  of  Engineers  (e.g.,  EM  2703,  EM 
2101)  and  nationally  recognized  specifications  and  standards  (e.g.,  ASCE 
Standard  7-88,  1990;  AISC,  1986)  specify  appropriate  nominal  load  values  and 
combinations  of  loads  to  be  used  in  design.  These  loading  requirements  have 
evolved  gradually  and  some  have  been  in  use  for  decades.  The  fact  that 
navigation  structures  generally  have  performed  well  indicates  that  the  current 
load  combination  requirements  are  conservative,  although  to  an  unknown  degree. 

Events  giving  rise  to  structural  loads  are  uncertain  in  their  magnitude  and 
structural  effect,  and  the  loads  vary  randomly  in  space  and  in  time.  These 
uncertainties  arise  from  inherent  randomness  in  naturally  occurring  or 
man-made  phenomena  and  from  assumptions  made  in  modeling  and  predicting 
behavior  of  complex  engineered  facilities.  Let 

U(t)  =  Xi  (t)  -h  X2  (t)  +  ...  -hXn  (t)  (1.1) 

represent  a  combination  of  structural  load  effects  (or  structural  actions  due  to 
randomly  occurring  operational  and  environmental  events),  as  illustrated  in 
Figure  1.1.  The  maximum  of  the  combined  load  effects  in  time  interval  o<t<T  (in 
which  T  might  be  one  year,  the  service  life  of  the  facility,  or  its 
inspection/maintenance  interval,  depending  on  the  purpose  of  the  event 
combination  analysis)  is  denoted, 

Umax  =  max  U(t)  (1.2) 

o<t<T 

Prudent  design  would  require  that  the  facility  be  designed  for  some  appropriate 
and  conservative  fractile,  Ud,  of  the  cumulative  probability  distribution 
function  (cdf)  of  Umax-  In  other  words,  we  would  like  to  assert  that  the 
probability  of  exceeding  the  design  envelope  furnished  by  Ud  during  the  service 
life  of  the  facility  is  acceptably  small. 

The  extreme  values  of  individual  load  effects  rarely  occur  simultaneously 
unless  the  peak  loads  are  associated  with  the  same  event  (see  Figure  1.1). 
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Engin0©rs  and  d0cision-mak6rs  g0n0rally  hav©  r©cogniz©d  intuitiv0ly  that 
combining  th©  p©ak  load  ©ff©cts  for  d©sign  purpos©s  is  unduly  cons©rvativ©. 

This  r©cognition  pr©dat©s  any  rational  consid©ration  of  ©v©nt  combination 
analysis  by  probabilistic  methods.  In  general, 

Umax  <  (t)  +  max  X2  (t)  +  .  .  . 

Traditional  practice  has  been  to  apply  "load  combination  factors"  less  than  unity 
in  estimating  the  design  value  of  the  combined  load,  Ud-  For  example,  Section 
2.3  of  ASCE  7-88  (ASCE  1990)  for  allowable  stress  design  requires  that  for 
combinations  of  live  and  wind  load, 

Ud  =  0.75(D -h  L -H  W) 

in  which  D,  L  and  W  =  dead,  live  and  wind  load  from  ASCE  7-88.  The  ASD 
Specification  for  steel  structures  (AISC,  1989)  addresses  this  problem  in  a 
slightly  different  way,  increasing  the  allowable  stress  by  a  factor  of  1/0.75  or 
1.33  for  combinations  involving  wind  or  earthquake  load  but  leaving  the  load 
combination  as  D  +  L  +  W. 

There  is  an  unstated  small  but  finite  (and  presumably  acceptable) 
probability  in  all  th©  above  cases  that  the  load  or  combined  load  effect  used  for 
design  will  be  exceeded  in  a  given  year.  The  fact  that  each  load  or  load 
combination  has  a  small  probability  of  being  exceeded  and  that  peak  values  of 
loads  in  combination  do  not  combine  concurrently  has  been  recognized  by  many 
standard  writers,  but  often  is  not  stated  explicitly  in  the  building  code  or 
commentaries.  Instead,  the  uncertainties  that  arise  from  inherent  variability 
and  imperfect  modeling  are  taken  into  account  through  the  us©  of  judgmental 
factors  of  safety.  Usually  there  is  no  need  for  the  engineer/designer  to 
consider  such  probabilities  explicitly  as  part  of  the  design  process.  Indeed,  the 
regulatory  climate  in  the  United  States  would  make  it  difficult  to  implement  a 
fully  probabilistic  design  procedure  in  which  designers  choose  the  probability 
level  depending  on  the  performance  requirements  for  the  facility. 

During  the  past  two  decades,  however,  substantial  advances  have  been 
made  in  probability-based  event  and  load  modeling  and  in  acquiring  the 
necessary  supporting  data  to  implement  probability-based  criteria  for  many 
phases  of  engineered  facility  design.  It  now  is  possible  in  many  instances  to 
^0't@i’fnin0  appropriate  design  loads  and  load  combinations  that  are  based  on  a 
specified  probability  of  being  exceeded  (Ellingwood,  et  al,  1982,  Galambos,  et 
al,  1982).  For  example,  the  nominal  loads  in  ASCE  7  -88  are  specified  as  N-year 
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mean  recurrence  interval  (MRI)  values,  or  the  value  corresponding  to  a 
probability  of  1/N  of  being  exceeded  in  a  given  year.  Parameter  N  customarily 
is  taken  as  50  years  for  wind  or  snow  loads  acting  on  ordinary  buildings;  a 
50-year  MRI  windspeed,  for  example,  has  a  2%  probability  of  being  exceeded 
annually.  The  earthquake  zone  map  is  based  on  a  seismic  hazard  analysis  which 
determines  contours  of  ground  acceleration  with  a  10%  probability  of  being 
exceeded  in  50  years;  the  corresponding  MRI  is  475  years  and  the  annual 
probability  is  0.21%  (FEMA,  1992). 

Structural  reliability  theory  and  probabilistic  load  combination  analysis 
methods  also  lead  to  better  insight  and  models  of  how  events  and  structural 
loads  may  actually  combine.  These  methods  have  been  applied  successfully  to 
building  structures,  reflected  by  the  proposed  design  load  combinations  in  ASCE 
7-88  for  load  and  resistance  factor  design  for  combinations  involving  live. 


snow,  wind  and  earthquake  effects: 

Ud  =  1.2  D  -H  1.6  L  (1.5) 

Ud  =  1.2  D  -h  1.6  S -f- 0.5L  +  0.8W  (1.6) 

Ud  =  1.2  D  -H  1 .3  W -I- 0.5  L  (1.7) 

Ud  =  1.2  D  +^.QE  +  0.5  L-h  0.2  S  (1.8) 

Ud  =  0.9  D  -H  (1.3W  or  1.0  E)  (1.9) 


(The  earthquake  load  combinations  in  Eqns  1.8  and  1.9  have  changed  recently  as 
part  of  the  move  to  implement  the  NEHRP  Recommended  Provisions  for  Seismic 
Regulations  (FEMA  1992).  Equation  1.9  is  used  to  check  overall  stability  of  the 
structure  or  situations  where  the  dead  and  wind  or  earthquake  effects 
counteract  one  another.  The  load  factors  0.5  on  L  or  0.2  on  S  in  combination 
with  other  effects  reflects  the  small  probability  that  extreme  values  of  the 
two  loads  occur  simultaneously. 

On  the  material  side,  the  AISC  LRFD  Specification  also  is  based  on 
probabilistic  reasoning.  For  example,  the  design  strength  for  compact  steel 
members  subjected  to  axial  tension  (yield),  axial  compression  and  bending  are 
given  by  (AISC  1986), 

0.90FyAg  >  Ud  (1.10) 
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0.85F  crAg 

>  Ud 

(1.11) 

0.90FyZx 

>  Ud 

(1.12) 

in  which  Fy  =  yield  strength,  Fcr  =  critical  strength,  Ag  =  gross  section  area,  and 
Zx  =  plastic  section  modulus,  and  Ud  represents  the  structural  action  that  is 
dimensionally  consistent  with  the  design  strength.  In  sum,  the  reliability  basis 
for  engineering  decision  analysis  allows  consistency  in  reliability  and 
performance  to  be  built  into  the  design  requirements  that  is  not  possible  when 
design  loads  and  safety  factors  are  selected  judgmentally. 

The  development  of  probability-based  event  combination  analysis 
guidelines  for  steel  or  reinforced  concrete  navigation  structures  operated  and 
maintained  by  the  U.S.  Army  Corps  of  Engineers  has  only  recently  been 
attempted,  with  an  application  to  steel  miter  gates  for  lock  and  dam  structures 
(Ellingwood  1993).  Earthquakes,  floods  and  similar  environmental  events  with 
the  potential  to  damage  a  structure  severely  are  rare  events  in  comparison  with 
events  giving  rise  to  design-basis  operating  loads.  Since  the  occurrences  of 
such  events  generally  are  statistically  independent,  the  probability  of  a 
coincidence  of  two  (or  more)  of  them  is  very  low  indeed.  In  such  cases,  any 
requirement  that  such  events  be  combined  in  design  may  contribute  little  to 
long-term  reli^Iity  and  perfdTmance^df  navigatiorT  aM  as  a  result  may  impose 
an  severe  economic  burden  on  the  Corps  and  the  public  with  little  benefit  to 
public  safety.  Progress  in  probability-based  event  combination  analysis  has 
been  made  in  other  contexts,  e.g.,  nuclear  power  plant  systems  (Shinozuka,  et  al, 
1981;  Schwartz,  et  al,  1981),  offshore  platforms,  and  other  critical  facilities. 
The  application  of  this  methodology  to  navigation  structures  remains  to  be 
explored. 

1.2  Research  Objectives  and  Scope 

This  report  develops  an  improved  methodology  for  determining 
design-basis  events,  loads  and  load  combinations  for  use  in  design  of  navigation 
structures  operated  by  the  U.S.  Army  Corps  of  Engineers.  This  methodology 
provides  a  rationale  to  support  the  implementation  of  design  guidelines  for 
event  combinations  in  which  the  events  are  taken  at  less  than  their  full  design 
intensity,  e.g.,  combinations  of  operating  and  flood  loads.  This  concept  also  has 
been  proposed  for  implementation  in  the  model  Eurocodes  (1991). 
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2.  PROBABILISTIC  MODELS  OF  LOAD  PROCESSES 


Probabilistic  event  or  load  combination  analysis  is  data-intensive  and 
requires  information  to  describe  the  rate  of  occurrence,  duration,  and  intensity 
of  all  loads  affecting  facility  performance.  This  section  identifies  basic 
mathematical  models  of  events  and  structural  loads  and  describes  the 
collection  and  analysis  of  available  data  to  identify  those  temporal 
characteristics  and  intensities  that  are  pertinent  to  the  study  of  their 
combinations  for  designing  navigation  structures. 

2.1  Basic  Load  Process  Models 

Probabilistic  descriptions  of  operating  and  environmental  events  and 
structural  actions  (forces,  deformations)  arising  from  such  events  must  convey, 
as  a  minimum,  a  description  of  the  rate  of  occurrence  and  duration  of  the  event 
and  the  intensity  of  the  structural  action.  The  occurrence  of  many  events  giving 
rise  to  time-dependent  structural  loads  can  be  modeled  as  a  Poisson  process 
(Larrabee  and  Cornell,  1981;  Pearce  and  Wen,  1984).  With  this  model,  the 
occurrence  of  loads  or  change  in  load  magnitudes  is  assumed  to  be  described  by 
the  Poisson  probability  law  (Papoulis,  1965); 

P[N(t)  =  n]  =  (X  t)n  exp(-X  t)/n!  ;  n  =  0,1,2,...  (2.1) 

in  which  N(t)  =  number  of  events  to  occur  during  interval  (0,t).  Three  sample 
functions  of  a  Poisson  process  are  illustrated  in  Figure  2.1.  In  Figure  2.1(a),  the 
event  is  more  or  less  continuous  and  the  load  is  always  present  (Type  I).  In 
Figure  2.1(b),  the  event  and  the  load  are  intermittent  (Type  II).  In  Figure  2.1(c), 
the  duration  of  the  event  is  neglected,  leading  to  a  sequence  of  load 
impulses. Each  load  event  is  characterized  by  a  mean  rate  of  occurrence,  X,  a 
mean  duration,  x,  and  a  cumulative  distribution  function  (cdf),  F(x),  modeling 
pulse  intensity.  Assuming  that  the  successive  pulse  intensities  are  identically 
distributed  and  statistically  independent  random  variables,  a  knowledge  of  the 
mean  rate  of  occurrence  of  the  events,  their  duration,  and  the  cdf  of  the  pulse 
intensity  is  sufficient  to  characterize  the  stochastic  load  process. 

If  the  pulse  process  is  always  nonzero  (Type  I),  the  parameters  X  and  x  are 
related  by  X  =  1/x  and  the  duration,  Td,  of  each  event  is  an  exponential  random 
variable  with  expected  value,  E[Td]  =  x.  If  the  pulse  process  is  intermittent 
(Type  II),  the  probability,  p,  that  the  process  is  nonzero  at  any  time  can  be 
obtained  from. 
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time 


( a )  Sustained  event  or  load  ( Type  I ) 


(2.2) 


N(T) 

P=  ^STdk=  ^  (:^T)  E[Td]  = 

oo  '  k  T-^  oo  ^ 

The  Poisson  load  process  model  has  been  shown  in  numerous  studies  to 
provide  an  adequate  representation  of  many  common  structural  loads  when  the 
structural  action  during  the  load  event  has  essentially  a  constant  or  static 
effect.  When  the  structural  action  is  dynamic,  the  occurrence  of  the  load  event 
still  can  be  described  by  a  Poisson  process;  however,  the  structural  load  may 
fluctuate  during  the  event  occurrence,  and  other  techniques  are  required  to 
analyze  its  structural  effect.  Fluctuating  loads  due  to  wind  or  earthquake 
sometimes  can  be  approximated  conservatively  as  Poisson  pulse  processes  if 
one  imagines  that  the  intensity  of  the  pulse  corresponds  to  the  maximum 
structural  action  that  occurs  during  the  event.  When  time-dependence  must  be 
considered,  (the  extremes  of  the  fluctuations  can  be  described  by  a  mean 
upcrossing  analysis  (Veneziano,  et  al,  1977),  as  described  subsequently. 

Time-varying  responses  X(t)  arising  from  randomly  occurring  events  can 
be  characterized  in  a  first-order  sense  by  their  pdf  Fx(x)  or  cdf  fx(x)  and  a  mean 
upcrossing  function,  Vx(x),  defined  as  the  rate  at  which  X(t)  crosses  level  x 
from  below.  This  is  illustrated  in  Figure  2.2.  In  the  special  cases  when  the 
changes  in  the  load  intensity  are  described  by  a  Poisson  process  with  rate 
parameter  A,,  and  the  load  intensities  are  described  by  a  sequence  of  identically 
distributed  and  statistically  independent  random  variables  with  cdf  Fs(x)  and 
pdf  fs(x),  Vx(x)  can  be  obtained  in  closed  form  (Winterstein  and  Cornell,  1984). 
For  example,  when  X(t)  is  always  nonzero  (Figure  2.1a): 

fx(x)  =  fs(x) 
vx(x)  =  X  Fs(x)Gs(x) 

in  which  Gs(x)  =  1-Fs(x)  =  complementary  cdf.  For  the  intermittent  pulse 
process  (Figure  2.1b): 

fx(x)  =  Xt  fs(x)  +  (1  -  At)  5(x)  (2.4a) 

vx(x)  =  A  Gs(x)  [1  -  At  Gs(x)]  (2.4b) 


(2.3a) 

(2.3b) 
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in  which  5(x)  =  Dirac  delta  function.  Finally  for  the  impulse  process  (Figure 
2.1c): 


vx(x)  =  X  Gs(x)  (2.5) 

More  generally,  the  upcrossing  rate  function  can  be  evaluated  from  the 
relation 


Vu(x)  =  /P[x(t)  <x  and  x(t+At)>x]\ 


(2.6) 


If  the  load  process  is  sufficiently  well-behaved  that  vx(x)  can  be  estimated 
from  one  sample  of  data  (mathematically,  X(t)  is  ergodic),  and  if  a  reasonably 
long  sample  of  data  is  available,  then  vx(x)  can  be  estimated  directly  from  this 

sample  of  data  using  Eqn  2.6  directly.  If  X(t)  and  its  derivative  X(t)  are 
continuous,  the  mean  upcrossing  rate  can  be  expressed  from  the  joint  density 
function  of  X(t)  and  X(t),  defined  asfx,x(x,x),  as  (Rice,  1945), 


vx  (x)  = 


X  fx  x(x,x)  dx 


(2.7) 


provided  that  the  joint  density  function  is  available.  Unless  X(t)  is  Gaussian, 
however,  this  joint  density  function  is  unknown  or  difficult  to  ascertain. 


In  addition  to  its  use  in  describing  local  fluctuations  in  load  intensity,  the 
mean  upcrossing  rate  function  plays  a  central  role  in  determining  the 
probability  that  X(t)  exceeds  a  barrier  level  x  during  the  interval  of  time  (0,  T), 
as  will  be  shown  subsequently. 

The  load  statistics  in  the  following  sections,  particularly  those  related  to 
environmental  load  effects,  are  presented  at  a  level  of  detail  necessary  to 
perform  the  event  combination  analysis  in  subsequent  sections.  A  summary  is 
presented  in  Table  2.1.  Environmental  load  statistics  are  strongly 
site-dependent.  Additional  published  data  are  available  for  those  wishing  to 
perform  a  site-specific  analysis,  as  noted  in  the  references. 
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2.2  Permanent  Loads 


2.2.1  Dead  Loads 

Dead  loads  are  random  in  magnitude  but  are  fixed  in  position  and  usually 
invariant  in  time  if  remodeling  effects  are  neglected.  They  arise  from  the 
weight  of  the  structure  and  permanent  equipment  and  attachments.  The  dead 
load  can  be  assumed  to  be  described  by  a  normal  cumulative  probability 
distribution  function  (c.d.f.),  with  a  mean  value  of  1.0  to  1.05  times  the  nominal, 
D,  and  a  coefficient  of  variation  (c.o.v.)  equal  to  0.10.  A  major  source  of  this 
variability  arises  from  the  treatment  of  point  loads  as  area-averaged  loads  for 
design  purposes. 

2.2.2  Temperature  and  other  Self-Straining  Effects 

In  ordinary  building  structures,  temperature  influences,  shrinkage  and 
creep  are  generally  considered  only  when  they  may  result  in  forces  or 
deformations  that  may  be  significant  to  the  ultimate  limit  state  considered.  In 
the  upper  Mississippi  region,  the  maximum  average  daily  temperature  variation 
{during  the  summer  months)  is  about  25°  F.  It  is  believed  that  forces  from  this 
magnitude  of  variation  are  insufficient  to  be  a  significant  design  consideration. 
Shrinkage  and  creep  can  be  controlled  by  proper  detailing  of  reinforcement. 

2.2.3  Prestressing 

Prestressing  is  a  type  of  imposed  deformation.  Current  codes  require  that 
prestressing  tendons  be  stressed  to  0.7  fpu  immediately  following  prestress 
transfer.  There  are  only  limited  statistical  data  on  prestressing  forces.  The 
cdf  of  prestress  effect  may  be  assumed  to  be  normal,  with  a  mean  that 
generally  is  close  to  the  nominal  value,  and  a  c.o.v.  of  about  0.02  at  transfer  and 
about  0.05  after  all  losses  have  occurred  (Ellingwood,  1983).  In  view  of  the 
magnitude  of  these  statistics  in  comparison  with  others,  it  is  reasonable  to 
assume  that  prestressing  effects  are  deterministic  in  concrete  structures  at 
navigation  facilities. 

2.3  Environmental  Loads 
2.3.1  Snow  Loads 

Snow  loads  are  significant  for  the  design  of  roofs  in  the  northern  tier  of 
states.  Stochastic  models  of  snow  load  on  roofs  are  developed  using 


meteorological  data,  supplemented  by  limited  snow  surveys  (Ellingwood  and 
Redfield,  1983;  O'Rourke,  et  al,  1982).  Snow  loads  are  dependent  on  site  and 
roof  geometry  exposure  (ASCE  7  1990)  and  are  seasonal  in  nature,  in  contrast 
with  other  common  environmental  loads. 

The  benchmark  statistics  needed  to  construct  the  snow  load  pulse  process 
are  derived  from  the  annual  extreme  snow  load  because  this  is  where  the  focus 
has  been  in  the  data  collection  and  analysis  that  supports  standard  development 
(ASCE  7  1990).  The  annual  extreme  roof  snow  load  can  be  described  by  a 
lognormal  distribution,  with  a  mean  value  that  typically  is  about  0.2S  and 
a  c.o.v  that  typically  is  about  0.90.  The  mean  rates  of  occurrence  of  significant 
snows  during  the  winter  season  are  site-dependent.  In  temperate  climates, 
there  may  be  periods  of  melting,  and  the  load  process  is  intermittent  (Type  II  in 
Figure  2.1(b)).  The  mean  rate  of  occurrence,  vs,  during  the  snow  season  would 
be  on  the  order  of  2  -  4  pulses/yr  with  a  mean  duration  of  perhaps  1  week.  In 
northern  climates  or  at  high  elevations,  the  snow  accumulates  during  the  winter 
season,  and  the  snow  load  can  be  modeled  (conservatively)  with  1  pulse/yr  and  a 

duration  which  may  be  as  much  as  3  months.  The  cdf  of  the  annual  extreme 

snow  load,  Fs(x),  is  related  to  the  cdf  of  the  pulse  intensity,  Fi(x),  by, 

Fs(x)  =  exp  -  F(x)]]  (2.8) 

Thus,  the  cdf  of  the  pulse  intensity  can  be  recovered  from  a  knowledge  of  the 
annual  extreme  and  temporal  characteristics  by, 

F(x)  =  1  -h  (1/Xs)  In  Fs(x)  (2.9) 

At  the  upper  fractiles  of  the  cdf  that  are  of  interest  in  structural  reliability 
analysis,  Fs  and  F  can  be  modeled  by  Type  I  distributions  of  largest  values 
(Benjamin  and  Cornell,  1970): 

F(x)  =  exp  [  -  exp  [  -  a(x  -  u)]  ]  (2.10) 

in  which  a,  u  =  parameters  of  the  distribution. 

2.3.2  Wind  Load 

Wind  loads  depend  on  the  wind  speed  and  building  exposure  and  geometry 
and  are  intermittent  in  nature  (Type  II  in  Figure  2.1b).  The  wind  speeds  used  for 
design  are  determined  from  basic  climatological  data,  while  the  building 


aerodynamic  effects  are  taken  into  account  by  wind  tunnel  modeling, 
supplemented  by  limited  in-situ  pressure  measurements. 

As  with  the  snow  load,  the  benchmark  statistics  used  to  model  the  wind 
load  as  a  pulse  process  are  derived  from  the  annual  extreme  wind  load  because 
this  is  where  the  focus  has  been  in  the  data  collection  and  analysis  that  support 
standard  development  (ASCE  7  1990).  The  statistical  data  are  site-dependent. 
The  annual  extreme  wind  load  can  be  modeled  with  a  Type  I  distribution  of 
largest  values,  with  a  mean  value  typically  about  0.30W,  in  which  W  =  nominal 
wind  load  (based  on  a  50-yr  mean  recurrence  interval  wind  speed)  and  a  c.o.v 
typically  about  0.60  (Galambos,  et  al,  1982).  The  mean  rate  of  occurrence  of 
structurally  significant  winds  is  about  2  -  4/yr  in  extratropical  regions,  with 
an  average  storm  duration  of  approximately  4  hours.  The  cdf  of  the  pulse 
intensity  (maximum  fluctuating  wind  load  during  the  occurrence  of  one  storm) 
is  defined  by  Eqn.  2.9,  substituting  wind  load  parameters  for  those  of  snow  load. 

2.3.3  Earthquake  Effects 

There  has  been  substantial  research  in  earthquake-resistant  design  since 
the  San  Fernando  earthquake  of  1971.  Most  of  this  research  has  shown  that  the 
earthquake  hazard  tends  to  be  underestimated,  particularly  in  the  Eastern 
United  States.  The  recurrence  intervals  between  large  events  in  the  Eastern  US 
are  large  and  as  a  result  there  is  little  or  no  experience  with  earthquakes  and 
their  social  impact  in  that  part  of  the  country.  Nontheless,  their  damage 
potential  should  be  considered  in  designing  navigation  structures  because  of 
their  consequences  for  the  economy  and  for  public  safety. 

Earthquake  loads  should  be  modeled  as  intermittent  (Type  II  in  Figure 
2.1b)  in  nature.  Earthquake  effects  on  building  structures  usually  are 
determined  using  an  equivalent  static-for-dynamic  analysis  in  which  the  base 
shear  is  computed  from  a  ground  motion  parameter  (usually  acceleration,  but 
sometimes  velocity)  and  general  dynamic  structural  characteristics.  The 
lateral  forces  that  are  statically  equivalent  to  this  base  shear  are  distributed 
over  the  structure  in  a  manner  consistent  with  its  mode  shape.  In  the  new 
NEHRP  provisions  (FEMA  1992),  the  (modal)  base  shear  is. 


V  =  CsW 


(2.11) 


in  which  Cs  =  seismic  design  coefficient  and  W  =  represents  the  effective 
mass-related  weight  of  the  structure.  Coefficient  Cg  is  defined  as, 

Cs  =  1.2AvS/RT2/3 


(2.12) 


in  which  Ag  and  Av  =  effective  peak  acceleration  and  velocity-related 
accelerations  determined  from  seismic  ground  acceleration  maps  or 
site-specific  seismic  hazard  analysis,  S  =  coefficient  for  soil  profile 
characteristics,  R  =  structural  response  modification  factor,  and  T  = 
fundamental  (modal)  period  of  the  structure.  The  coefficient  Cg  is  tantamount 
to  an  inelastic  yield  spectrum  for  an  oscillator  with  5%  damping. 

For  hydraulic  structures,  the  inertial  hydrodynamic  forces  are  determined 
from  Westergaard's  equation  (Westergaard,  1933), 

p(y)  =  0.875  Yw  Av  Hy  (2.13) 

in  which  yw  =  density  of  water,  Hy  =  water  depth,  and  y  =  depth  below  water 
surface. 

The  value  of  ground  acceleration  or  velocity  (or  spectral  acceleration) 
used  in  design  is  determined  from  a  seismic  hazard  analysis.  Seismic  hazard 
analysis  carries  with  it  the  largest  component  of  uncertainty  in  the  process 
leading  to  earthquake  forces  on  engineered  facilities  (Reiter,  1990).  Research 
in  earthquake-resistant  design  during  the  past  15  years  has  enabled  seismic 
hazard  analysis  and  design  ground  motions  at  a  particular  site  to  be  placed  on  a 
probabilistic  basis  (Algermissen,  et  al  1982).  Earthquake  hazards  in  the 
Western  United  States  generally  can  be  associated  with  a  series  of  capable 
faults.  Such  faults  are  not  apparent  in  the  Eastern  United  States,  and  the 
seismic  hazard  analysis  there  begins  with  an  identification  of  postulated 
seismic  source  zones.  Seismicity  within  these  zones  of  potential  future 
earthquake  occurrence  is  determined  and  mean  rates  of  occurrence  of 
earthquakes  of  various  magnitudes  (or  intensities)  are  identified.  Attenuation 
functions  relating  effective  peak  ground  acceleration  to  magnitude  and 
epicentral  distance  are  used  to  propagate  the  earthquake  ground  motion  from 
the  source  to  the  building  site.  Finally,  a  probability  distribution  of  effective 
peak  ground  accelerations  at  the  site  is  determined  by  summing  (integrating) 
over  ail  possible  earthquake  sources  and  magnitudes  consistent  with  each 
underlying  source  hypothesis.  The  result  is  usually  presented  as  a 
complementary  cdf  or  seismic  hazard  curve,  GA(a),  showing  the  annual 
probability  of  exceeding  a  specified  ground  acceleration,  a. 

Typical  seismic  hazard  curves  for  building  sites  in  the  Western  US  and  the 
Eastern  US  are  compared  in  Figure  2.3.  The  seismic  hazard  curves  can  be 
described  by  a  Type  II  distribution  of  largest  values  (Cornell,  1968), 


GA(a)  =  1  *  FA(a)  =  1  -  exp[-(a/u)Oc] 


(2.14) 


in  which  a  =  ground  acceleration  and  u  and  a  are  parameters  of  the  distribution. 
Parameter  a  determines  the  slope  of  the  hazard  curve.  The  curve  in  the  Eastern 
US  is  much  flatter  (smaller  a)  than  that  in  the  Western  US  (see  Figure  2.3);  this 
is  because  of  the  relatively  larger  uncertainty  associated  with  Eastern 
seismicity  due  to  the  absence  of  large  historical  events  during  the  period  of 
modern  instrumentation.  The  Western  US  hazard  curve  saturates  at  acceleration 
levels  on  the  order  of  1.2  to  1.8  times  the  design  earthquake;  in  the  Eastern  US, 
the  saturation  level  is  unknown,  but  is  believed  to  be  on  the  order  of  2  to  3 
times  the  design  earthquake. 

Parameter  a  in  Eqn  2.14  is  related  to  the  c.o.v.  in  annual  maximum 
effective  peak  acceleration.  Underlying  the  hazard  analysis,  on  which  the 
seismic  risk  map  found  in  Section  9  of  ASCE  7-88  is  an  overall  value  ot  =  2.3  for 
the  continental  United  States  (Galambos,  et  al,  1982);  the  associated  c.o.v.  is 
138%.  Later  hazard  analyses  (Algermissen,  et  al,  1982)  on  which  the  NEHRP 
Recommended  Provisions  for  Seismic  Regulations  (FEMA  1992)  are  based  have 
not  led  to  significant  reductions  in  this  enormous  c.o.v.,  although  some  variation 
from  site  to  site  is  evident  in  the  later  study.  The  factor  a  tends  to  be  larger 
for  sites  in  the  Western  US,  decreasing  from  about  5.5  (c.o.v.  =  0.28)  at  San 
Francisco,  CA  to  approximately  2.3  (c.o.v.  =  1.38)  at  Boston,  MA  and  Memphis,  TN. 

The  NEHRP  Recommended  Provisions  are  being  considered,  with  minor 
variants,  for  adoption  by  ASCE  Standard  Committee  7  for  its  revision  to  the 
Standard  on  Structural  Loads  and  by  the  Model  Codes.  The  design  earthquake 
envisioned  for  buildings,  as  reflected  in  the  design  ground  motion  contour  maps 
is  based  on  a  specified  probability  of  0.10  or  less  of  being  exceeded  in  50  years; 
this  corresponds  to  a  mean  recurrence  interval  of  about  475  years.  Earthquakes 
with  magnitudes  less  than  Ml  =  4.0  or  intensities  less  than  V  were  not 
considered  in  evaluating  seismicity  for  this  map  because  it  was  believed  that 
their  capability  for  causing  structural  damage  was  negligible.  This  probability 
is  much  less  than  the  annual  probabilities  of  exceeding  the  design  snow  or  wind 
loads  in  ASCE  7-88,  which  normally  are  0.02.  The  uncertainty  in  the  earthquake 
effect  is  vested  in  this  conservative  specification  of  the  structural  action,  E, 
due  to  the  earthquake.  As  a  result,  the  load  factor  on  E  is  set  equal  to  1.0  in  the 
NEHRP  Recommended  Provisions  and  in  the  revision  to  ASCE  Standard  7  (cf  Eqns. 
1.8  and  1.9)  rather  than  a  greater  value. 

The  distributions  of  the  annual  extreme  effective  peak  ground 
acceleration  (EPA)  and  the  pulse  intensity  can  be  related  using  Eqn  2.9, 
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Figure  2.3  Typical  seismic  hazard  curves 


substituting  the  earthquake  statistics  for  those  of  snow.  Thus,  the  maximum 
EPA  of  an  individual  earthquake  can  be  described  by, 

F(a)  =  1  +  (1/>.e)  In  FA(a)  »  1 -GA(a)/>bE  (2-15) 

If  ao  denotes  the  minimum  level  of  EPA  necessary  for  a  structurally  damaging 
earthquake,  F(ao)  =  0  and  Xe  =  -In  FA(ao).  Parameter  ao  typically  would  be  about 
0.05g.  The  parameter  Xe  varies  from  about  0.01 /yr  for  sites  in  the  Eastern  US  to 
about  0  10/yr  for  sites  in  the  Western  US,  although  a  few  sites  in  California 
indicate  values  of  ;^e  as  high  as  0.5.  The  mean  duration  of  ground  shaking  would 
be  in  the  range  of  10  to  60  sec,  depending  on  the  magnitude  of  the  earthquake. 

The  annual  extreme  earthquake  structural  action  can  be  modeled  by  a  Type 
II  distribution  of  largest  values  (Galambos,  et  al,  1982).  The  coefficient  of 
variation  in  the  load  effect  typically  is  0.8  -  1.3,  practically  all  of  which  is 
attributable  to  the  uncertainty  in  the  basic  seismic  hazard  at  a  building  site. 

2.3.4  Flood  Loads 

A  design  flood  carries  with  it  considerations  for  both  hydrological  risk 
and  facility  risk.  However,  facility  risk  depends  on  the  operation 
characteristics  of  the  facility,  and  thus  would  vary  for  different  navigation 

structures. 

Flood  risk  estimates  are  region-dependent.  Risk  assessment  of  facilities 
such  as  dams  has  been  part  of  water  resources  planning  and  management  for 
many  years  (Cochrane,  et  al.  1987).  Design  flood  levels  have  been  most 
commonly  expressed  in  terms  of  either  annual  probability  of  being  exceeded  or 
as  a  probable  maximum  flood  (PMF).  The  recommended  probability  distribution 
for  annual  flood  flow  is  the  log-Pearson  Type  III  (Interagency,  1982).  The 
annual  probability  often  is  selected  equal  to  0.01,  corresponding  to  a  mean 
recurrence  interval  of  100  years.  Historical  records  of  flood  events  in  most 
areas  of  the  United  States  seldom  extend  back  more  than  50  years.  Thus, 
extrapolations  beyond  annual  probabilities  of  0.01  have  limited  validity  and  may 
be  subject  to  significant  sampling  error.  The  PMF  is  generated  from  a 
mathematical  model  of  local  climatology  and  involves  the  conversion  of 
historical  storm  and  precipitation  data,  which  generally  are  more  complete  than 
flood  data,  to  surface  runoff  and  flood  levels.  The  PMF  is  derived  from  the 
probable  maximum  precipitation,  estimates  of  which  are  available  from  the 
National  Weather  Service  (Committee.  1983).  These  mathematical  models 
introduce  uncertainties  of  their  own  into  the  determination  of  flood  level.  The 


annual  probability  associated  with  a  PMF  is  unknown,  but  likely  is  less  than 
10-4/yr,  It  should  be  emphasized  that  the  PMF  does  not  have  a  probabilistic 
basis  like  the  100-year  MRI  flood;  however,  studies  have  indicated  that  the 
ratio  of  24-hr  PMP  to  100-year  mean  recurrence  interval  precipitation  at  sites 
in  the  upper  Mississippi  basin  is  approximately  5  (Safety,  1985). 

The  flood  statistics  of  interest  for  event  combination  analysis  are 
primarily  the  mean  rate  of  occurrence  and  the  duration  of  the  flood.  It  seems 
unlikely  that  a  navigation  facility  would  be  designed  to  withstand  the  effects  of 
extreme  floods  directly.  For  example,  the  gates  at  a  lock  and  dam  might  simply 
be  opened  during  a  design-basis  flood  to  allow  passage  of  high  levels  of  water 
without  damage  to  engineered  structures. 

2.4  Operating  Loads 

The  normal  operations  of  navigation  structures  such  as  gates  at  lock  and 
dam  facilities  give  rise  to  several  types  of  structural  loads.  These  include 
hydrostatic  pressure,  temporal  head,  hydrodynamic  load,  equipment  loads,  and 
impact.  The  probabilistic  modeling  of  these  loads  are  dependent  on  the 
operating  characteristics  of  the  facility.  Little  statistical  data  are  available 
to  support  these  probabilistic  load  models.  Statistical  analysis  of  poolwater 
and  tailwater  elevations  at  three  lock  and  dam  sites  on  the  Mississippi  River 
has  provided  some  data  on  hydrostatic  forces,  while  the  analysis  of  a  survey  of 
lock  operators  on  several  major  river  systems  has  provided  some  preliminary 
data  on  the  other  operating  loads.  The  analysis  of  these  data  is  described  in 
detail  elsewhere  (Ellingwood  1993);  here,  only  a  summary  is  provided. 

2.4.1  Hydrostatic  Loads 

The  hydrostatic  load  Hg  on  a  gate  in  the  closed  position  is  determined  by 
the  upper  and  lower  pool  elevations.  These  elevations  are  site-dependent  and 
their  variation  is  determined  on  weather  conditions.  While  there  is  natural 
variation  in  pool  elevation  from  year  to  year,  the  variability  in  elevation  and  in 
differential  is  relatively  small  because  the  river  level  is  controlled  by  dams  for 
flood  control  upstream  and  downstream.  The  resultant  force  per  unit  width  of 
structure  is, 

F  =  Yw  A  (Hp  -h  Ht)/2  (2.16) 

in  which  A  =  Hp  -  Ht  and  Hp,  Ht  =  pool  and  tailwater  depths.  The  hydrostatic 
pressure  at  depth  z  below  the  upper  pool  elevation  is. 
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Yw  z;  z  <  A 


(2.17) 


P(z)  =  { 

Yw  A;  z  >  A 

The  maximum  annual  differential  does  not  necessarily  coincide  with  the 
maximum  pool. 

The  randomness  in  the  intensity  of  F  or  p(z)  is  determined  primarily  by  the 
relative  control  on  water  level  provided  by  upstream  and  downstream  dams.  For 
Lock  and  Dam  Nos.  24,  25  and  27  on  the  Mississippi  River,  the  coefficients  of 
variation  in  F  ranged  from  less  than  0.10  to  about  0.30.  The  mean  rate  of 
occurrence  and  mean  durations  of  these  hydrostatic  loads  is  determined  by  the 
rate  of  lockages;  for  heavily  travelled  rivers  such  as  the  Mississippi  or  Ohio, 
typical  values  might  be  X,  =  200/month  and  x  =  90  minutes. 

2.4.2  Equipment  and  Hydrodynamic  Loads 

Hydrodynamic  load  Hd  and  equipment  load  Q  arise  during  the  operation  of 
the  gates.  Under  normal  operating  conditions.  Q  and  Hd  represent  a  statically 
equilibrated  force  system.  However,  if  there  is  a  submerged  obstruction,  the 
bottom  of  the  gate  leaf  may  bind  during  operation  while  the  operating  load  is 
applied  at  the  top  of  the  leaf,  causing  the  gate  leaf  to  twist.  Uncertainties  in  Q 
and  Hd  may  arise  from  variations  in  pool  elevations,  from  surge  effects  caused 
by  movement  of  the  gates  or  overtravel  of  water  in  the  culvert  system  during 
filling.  Temporal  hydraulic  loads  occur  due  to  temporary  variation  in  water 
level  and  wave  action  caused  by  tow  movement  or  from  overfilling  or 
underfilling  the  local  chamber.  The  load  processes  for  equipment,  hydrodynamic 
and  temporal  head  were  determined  from  a  Delphi  of  lock  and  dam  operators 
{Ellingwood  1993). 

2.4.3  Impact 

Impact  loads  arise  from  accidents  in  which  vessels  collide  with 
navigation  structures.  Such  accidents  may  result  from  failure  of  steering  or 
towing  mechanisms,  wind  and  current,  or  pilot  error.  A  database  on  impact 
occurrence  is  maintained  by  each  USAGE  District  Office.  These  data  are 
summarized  for  lock  and  dam  structures  in  TR-REMR-HY-7  (Martin  and  Lipinski, 
1990).  A  review  of  these  data  indicated  that  the  average  rate  of  accidents 
involving  miter  gates  at  locks  on  the  Mississippi,  Illinois  and  Ohio  Rivers  is 
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approximately  0.06/month;  however,  not  all  of  these  accidents  required  that  the 
lock  be  shut  down  for  repair.  When  only  those  accidents  that  involved 
significant  damage  to  the  lock  (e.g.,  over  $50,000  in  damage)  are  considered,  the 
mean  rate  of  occurrence  decreases  to  approximately  0.016/month,  or  about  once 
every  five  years. 

The  force  due  to  barge  impact  depends  on  the  mass  of  the  barge  and  its 
velocity  at  impending  collision.  Forces  on  navigation  structures  due  to  barge 
impact  can  be  enormous;  in  a  recent  impact  test  at  LD  26  on  the  Mississippi 
River  (Chasten  1992),  the  force  measured  as  a  fully  loaded  9-barge  tow 
impacted  the  gate  at  0.94  ft/s  (0.64  mph)  was  over  600  kips.  Terminal 
velocities  of  this  order  are  not  inconsistent  with  the  limited  data  available 
from  other  sources  (Ellingwood,  1993). 

2.5  Summary 

A  summary  of  the  event  and  load  process  characteristics  used  in  the 
subsequent  event  and  load  combination  analysis  is  presented  in  Table  2.1.  These 
parameters  include  the  mean  rate  of  occurrence  and  duration,  X  and  x,  and  the 
point-in-time  probability  distributions  F(x). 
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Table  2.1  -  Statistics  of  Structural  Load  Events 


Load 

Occurrence 
Rate(/yr)  Duration 

Dead 

const. 

50  yr 

Snow 

4 

1  wk 

Wind 

2 

4  hr 

Earthq. 

0.02 

30  s 

Operations 

4,800 

2  m 

Hydrost. 

2,400 

90  m 

Temporal 

4,800 

5  m 

Impact 

0.0013 

15  s 

Intensity 

Mean 

c.o.v. 

cdf 

1.05D 

0.10 

Normal 

0.20S 

0.90 

Type  1 

0.30W 

0.60 

Type  1 

Site-dep 

0.85 

Type  II 

0.70 

0.50 

Type  1 

0.9Hs 

0.20 

Type  1 

0.5Ht 

0.80 

Type  1 

1.01m 

n.a. 

n.a. 
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3.0  PROBABILISTIC  BASES  FOR  EVENT  COMBINATION  ANALYSIS 


3.1  General  Considerations 

In  performing  a  probabilistic  event  combination  analysis  with  a  view 
toward  ensuring  structural  safety,  the  primary  interest  is  in  structural  actions 
(forces,  deformations)  and  their  combinations  arising  from  such  events. 

Before  performing  such  an  analysis,  it  is  important  to  identify  any 
possible  dependencies  among  individual  loads  and  to  screen  from  consideration 
other  combinations  that  are  not  possible.  This  preliminary  evaluation  should  be 
done  to  simplify  the  event  combination  analysis,  which  can  be  data-intensive 
and  relatively  complex  mathematically.  For  example,  it  would  be  reasonable  to 
assume  in  evaluating  a  lock  and  dam  facility  that  gate  operation  is  independent 
of  the  environment  (leaving  aside  those  periods  of  the  year  where  winter 
weather  may  prevent  navigation  entirely  and  obviate  the  need  to  operate  the 
structure).  Thus,  one  would  assume  that  hydrostatic  pressure  and  earthquake 
could  be  treated  as  statistically  independent,  and  that  analysis  of  combinations 
of  extreme  snow  (or  ice)  and  operating  loads  would  not  be  required.  Similarly, 
if  one  were  considering  loads  acting  on  a  moveable  gate,  hydrostatic  force  Hs 
is  exclusive  of  operating  loads  Hd  or  Q  because  the  gate  does  not  operate  under 
differential  head;  on  the  other  hand,  Hd  and  Q  form  a  statically  equivalent  couple 
during  normal  operation. 

Consider  now  those  events  or  structural  actions  that  may  need  to  be 
combined  according  to  the  above  preliminary  assessment.  Most  events  giving 
rise  to  significant  structural  loads  occur  infrequently  and  usually  do  not  last 
very  long.  The  appropriate  combination  of  structural  actions  will  depend  in  a 
general  sense  on  the  rate  of  occurrence  of  the  events,  their  duration,  the 
magnitude  of  the  structural  action  and,  in  the  case  of  a  dynamic  force,  its  scale 
of  fluctuation  during  the  event.  Mathematically,  such  events  can  be  modeled  in 
the  macro-timescale  by  an  intermittent  pulse  process  of  some  kind  as  described 
in  Section  2.1,  one  that  captures  the  randomness  in  occurrence  and  duration  of 
the  events  and  a  general  measure  of  their  intensity.  If  the  structural  action  is 
essentially  static  and  varies  little  within  a  given  event  occurrence,  such  a 
pulse  process  model  may  be  sufficient  since  all  that  is  required  for  safety 
evaluation  is  the  extreme  load  to  occur  during  the  event.  For  example,  loads  due 
to  hydrostatic  pressure  or  snow  can  be  modeled  by  pulse  processes.  On  the 
other  hand,  if  the  structural  action  during  an  event  is  dynamic,  (e.g.,  earthquake, 
surge),  the  macroscale  modeling  of  the  load  process  may  not  be  sufficient,  and 
additional  steps  are  required  to  model  the  variations  in  the  micro-timescale. 
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Methods  for  modeling  such  effects  for  individual  loads  were  presented  in 
Section  2.1. 

3.2  Fundamental  Event  Combination  Analysis 

The  basis  for  the  event  combination  methodology  is  the  probability 
distribution  of  the  maximum  of  a  sum  of  time-  dependent  structural  loads 
modeled  as  stochastic  processes.  Let 


U(t)  =  X,(t)  +  X2(t)  +  .  .  . 

(3.1) 

Umax  =  max  U(t) 

(3.2) 

o<t<T 


in  which  Xi(t)  =  structural  action  due  to  a  stochastic  event  (e.g.,  operating  load, 
temporal  hydrostatic  load,  earthquake). 

If  two  or  more  loads  that  are  always  nonzero  are  combined,  the 
probability  distribution  of  the  combined  load  involving  Xi(t),  X  2(t),.  .  .,  at  any 
point  in  time  can  be  obtained  by  the  convolution  theorem  of  probability  theory 
{Papoulis,  1965).  For  example,  if  U  =  Xi  +  X2  ,  the  cumulative  probability 
distribution  function,  Fu(x),  is. 


Fu(x)  =  [  Fxi(X  -  X2)  fx2(X2)  clx2 


(3.3) 


in  which  Fxi(x)  is  the  cdf  of  X  1  and  fx2  (x)  is  the  density  function  for  X2. 

Analyzing  combinations  of  intermittent  loads  which  may  be  zero  for 
significant  periods  of  time  is  more  complex.  The  possibility  must  be  considered 
that  at  any  given  time,  only  one  of  the  loads  acts  or  that  two  (or  more)  of  the 
loads  act  simultaneously.  Moreover,  one  generally  is  more  interested  in  the 
maximum  of  the  combined  effects  during  a  service  period,  T,  rather  than  the 
point-in-time  combined  load  effect  for  design  purposes.  The  analysis  of  the 
maximum  of  a  combination  of  time-dependent  stochastic  events  is  related  to 
the  difficult  first-passage  problem  of  probability  theory.  Although  the 
distribution  of  Umax  generally  cannot  be  determined  exactly,  conservative 
approximations  are  available  (Larrabee  and  Cornell,  1981;  Pearce  and  Wen, 

1984;  Winterstein  and  Carroll  1984). 
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It  is  assumed  in  the  following  that  the  structural  actions  Xi(t)  represent 
stresses,  strains,  deformations  or  internal  forces  arising  from  event  i.  The 
limit  state  equation  is  written  in  terms  of  the  processes,  Xi(t),  rather  than  in 
terms  of  the  events.  However,  the  temporal  nature  of  the  event  and  the  load 
generally  are  directly  related.  Thus,  the  processes  Xi(t)  can  be  classified  into 
the  same  basic  three  categories  as  shown  in  Figures  1.1  and  2.1.  We  seek  a 
conservative  approximation  to  the  cdf  describing  Umax- 

Upcrossing  analysis  provides  a  powerful  tool  for  analyzing  extremes  of 
combinations  of  structural  loads(Veneziano,  et  al,  1977;  Rackwitz,  1985).  The 
load  process  can  be  nonstationary  in  microscale  without  significant  error  as 
long  as  it  is  stationary  in  macroscale.  It  can  be  shown  (e.g.,  Shinozuka,  1981) 
that  a  conservative  upper  bound  to  the  complementary  cdf  of  Umax  is  given  by. 


P[Umax  >  x]  <  P[U(0)  >  X  ]  +  Vu(x,t)dt  (3.4) 

in  which  vu(x,t)  is  the  (time-dependent)  mean  rate  at  which  U(t)  upcrosses  level 
X.  If  the  process  U(t)  is  stationary,  the  mean  upcrossing  rate  function  is 
independent  of  t,  and 

P[Umax  >  x]  <  P[U(0)  >  x]  +  Vu  (x)  T  (3.5) 

If  the  upcrossings  of  level  x  are  rare  as  they  usually  are  in  structural 
reliability  analyses,  the  occurrence  of  upcrossings  can  be  modeled  as  a  Poisson 
process  with  intensity  vu(x).  In  this  case,  the  complementary  cdf  is, 

P(Umax  >  x]  =  1  -  exp[-vu(x)T]  (3.6) 

Eqn.  3.6  is  asymptotically  correct  as  x  becomes  large.  At  lower  values  of  x,  the 
upcrossing  events  are  not  Poisson;  however.  It  has  been  found  that  a  good 
approximation  at  moderate  levels  of  x  is  (Ditlevsen,  1981), 

P[Umax  >  x]  =  1  -  Fu(x)  exp[-vu(x)T/Fu(x)]  (3.7) 

Consider  for  illustration  a  case  where  there  are  two  time-dependent 
structural  actions,  so  that 

U(t)  =  Xi(t)  +  X2(t)  (3.8) 
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Assume  that  Xi(t)  and  X2(t)  are  statistically  independent.  There  are  four 
mutually  exclusive  and  collectively  exhaustive  ways  for  Xi  and  X2  to  occur: 


Eoo  =  {Xi(t)  =  0  and  X2{t)  =  0} 

E10  =  {Xi{t)  >  0  and  X2(t)  =  0} 

Eo2  =  {Xi(t)  =  0  and  X2(t)  >  0} 

Ei2  =  {Xi(t)  >  0  and  X2(t)  >  0} 

The  mean  upcrossing  rate  vu(x)  can 
expectation: 


(3.9a) 

(3.9b) 

(3.9c) 

(3.9d) 

be  determined  from  the  theorem  of  total 


vu  (x)  =  Z  vu|E,(x)  P[Ei]  (3.10) 

j 

in  which  i  =  10,  02,  and  12.  This  can  be  expressed  as, 

vu(x)  =  V10  (x)  P10  +  V02  (x)  P02  +  V12  (x)  P12  (3*1 1) 

in  which  vio,  pio  =  P[Eio]  are  the  mean  upcrossing  rate  and  probability  of 
occurrence  of  Xi  when  X2  =  0;  V12.  Pi 2  are  mean  upcrossing  rate  and  probability 
of  coincidence  of  Xi  +  X2,  etc.  Implementation  of  Eqn  3.11  requires  information 
on  the  probabilities  P[Ei]  and  mean  upcrossing  rates  Vi. 

The  probabilities  P[Ei]  depend  on  the  distribution  functions  of  the 
inter-arrival  times  of  Xi  and  X2  and  durations  of  these  events.  Under  the 
assumption  that  the  occurrence  of  individual  load  events  can  be  described  by  a 
Poisson  process,  with  rate  parameter  X.  the  inter-arrival  times, Ta,  for  the 
events  are  independent  random  variables  described  by  an  exponential 
distribution,  with  parameter  X‘ 


Fyg(x)  =  1  -  exp[  -  X\  x];  x>0  (3.12) 

The  expected  value  of  Ta,  E[Ta],  is  MX.  We  now  suppose  that  the  durations,  Td,  of 
the  event  also  are  exponential  with  parameter,  \i,  but  impose  the  additional 
condition  that  Td  cannot  be  longer  than  Tg.  The  duration  of  the  event  thus  is 
described  by  the  following  conditional  cdf: 
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^Td(x)  =  P[Td  <  x|Td  <  Ta] 


{3.13a) 


which  can  be  expressed  as, 


FTd(x)  =  1  -  exp[ -(>.+ ^)x];  x>0 


(3.13b) 


The  expected  value  of  Td  is  E[Td]  =  M{X  +m.).  The  probability  that  the  event  is 
acting,  or  "on",  is  p  =  E[Td]/E[Ta]  =  r/(1  +  r),  in  which  r  =  Xl\i,  while  the 
probability  that  the  event  is  not  acting,  or  is  "off,"  is  1  -  p  =  1/(1+r).  When  Xl\i 
1.0,  the  process  is  almost  always  "on";  as  1/|i-4  0,  one  obtains  a  process 
consisting  of  impulses  with  (essentially)  no  duration,  as  in  Figure  2.1c. 


When  two  events  with  the  above  temporal  characteristics  are  possible, 
there  will  be  times  at  which  there  are  none,  one,  or  two  events  acting  at  one 
time.  A  standard  result  from  renewal  theory  states  that  the  probabilities  of 
having  neither  Xi  nor  X2,  Xi  without  X2,  X2  without  Xi,  and  Xi  and  X2  together  (cf 
Eqns  3.9)  are, 

2 

Poo  =  "i/n  (1  +  fj) 

i>1 


P10  =  ri  Poo 
P02  =  f2  Poo 


(3.14) 


P12  =  rir2  Poo 

Extensions  to  more  than  two  events  are  straightforward  (Shinozuka,  1981).  The 
parameters  n  generally  will  be  much  less  than  unity  for  transient  events 
arising  from  infrequent  environmental  events  such  as  floods  or  earthquakes. 

The  mean  upcrossing  rates  of  Xi  and  X2  in  Eqn  3.11  can  be  computed 
analytically  (Eqns  2.3  -  2.5)  or  from  a  time  history  of  the  load  for  a 
well-behaved  process  (Eqn  2.6).  For  the  combined  event  Xi  +  X2,  the  point 
crossing  formula  (Larrabee  and  Cornell,  1981) 

Vi2(x)=j^  vi(x-y)  f2(y)  dy j^V2(x-y)  fi(y)  dy  (3.15) 

expressing  the  combined  upcrossing  rate  in  terms  of  the  marginal  pdfs  and 
upcrossing  rates  can  be  used. 
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An  alternate  procedure  for  event  combination  analysis  is  denoted  the  load 
coincidence  method  (Wen,  1977).  Using  the  load  coincidence  method,  the 
complementary  cdf  is  given  by, 

P[U„„>xl.  I1uTGi  +  II  XiiTGii  +  ...  (3.16) 

i  i  j 


in  which  X\  =  mean  rate  of  occurrence  of  the  ith  load  acting  alone,  Gj  = 
probability  that  load  i  acting  alone  exceeds  x,  X\i  =  mean  rate  of  occurrence  of  xi 
+  Xi  Gii  =  probability  that  Xi  +  Xj  exceeds  x,  etc.  Eqn  3.16  implies  that  the  mean 
upcrossing  rate  is  equal  to  the  occurrence  rate  of  the  combinations  multiplied 
by  the  probability  that  the  combination  exceeds  x,  given  that  it  occurs. 


The  similarities  of  Eqns  3.5  and  3.11  to  3.16  should  be  noted.  The 
equations  are  identical  if 


Vj{x)T  =  Gi(x) 


(3.17a) 


Pi  =  Xi 


{3.17b) 


These  identities  hold  if  the  processes  are  "sparse,"  i.e.,  X\x\  <  '\.  This  is 
characteristic  of  processes  describing  significant  operating,  environmental  and 
accidental  loads,  which  occur  infrequently  and  have  short  duration.  The 
accuracy  of  both  approaches  has  been  established  by  Monte  Carlo  simulation 
(Wen,  1977;  Wen,  1981;  Winterstein  and  Cornell,  1984).  Both  approaches  show 
good  agreement  at  the  levels  of  x  significant  for  structural  safety  studies. 

Figure  3.1  illustrates  the  analysis  of  two  intermittent  events  with  an 
event  tree  (Winterstein,  1980)  to  facilitate  better  visualization  of  the  event 
combination  methodology.  The  occurrence  of  the  extreme  load  events  is 
assumed  to  be  modeled  by  Poisson  processes,  and  thus  the  probability  of  load  1 
(or  load  2)  occurring  in  interval  of  time  (t,t+dt]  is  X,idt  (or  X2  dt),  i.e.. 


P[load  1  occurs  in  {t,t+dt]]  =  X^d\ 

This  probability,  in  itself,  does  not  guarantee  that  only  load  1  acts,  since  it  is 
possible  that  load  2  is  also  acting  at  the  time  load  1  acts.  To  ensure  that  the 
load  events  are  mutually  exclusive,  the  possibility  of  a  joint  occurrence  of 
loads  1  and  2  must  be  taken  into  account.  The  probability  that  load  1  occurs 


27 


Event  2  begins  during 
event  1  (coincidence) 


Event  1  occurs  in 
(t,t+At ) 


•Event  2  does  not  occur  during 
event  1  ( event  1  alone ) 


►No  occurrence  of  either  event 


Neither  event  occurs  in  (t,t+At) 


.Event  1  begins  during 
event  2  (coincidence) 


Event  2  occurs  in 
( t , t+At ) 


^ Event  1  does  not  occur  during 
event  2  ( event  2  alone ) 


Figure  3*1  Event  tree  for  analysis  of  event  combinations 
(after  Winterstein^  1980) 


and  load  2  does  not  occur  during  the  interval  in  which  Xi  acts  is  ?tidt(1-  A 

similar  analysis  must  be  done  for  load  X2  in  combination  with  load  Xi.  The 
probability  of  a  coincident  event  occurring  in  (t,t+dt]  is  ?iiX2dt(Ti  +  t2).  The 
mean  rates  of  occurrence  of  the  mutually  exclusive  events  are  then: 


^02  =  ■  ^12  =  ^2  "  ^12  ('Tl  +  '^2) 

X^2  =  h'^2  ('Cl  +  'C2) 


(3.18a) 

(3.18b) 

(3.18c) 


The  mean  durations  of  the  mutually  exclusive  events  are  Ti  and  X2  and,  for  the 
coincident  event, 


T12  = 'C1'C2/('C1  + 'C2)  (3.19) 

If  the  occurrence  of  one  event  influences  the  occurrence  of  a  second  event, 
the  two  load  processes  are  no  longer  statistically  independent  and  Eqns  3.18  and 
3.19  are  not  applicable.  The  Poisson  occurrence  model  on  which  the  event 
combination  analysis  is  based  can  be  modified  to  take  this  dependence  in 
occurrence  into  account. 

The  situation  can  be  visualized  by  the  load  process  sample  functions 
illustrated  in  Figure  3.2.  Suppose  that  event  1  occurs  according  to  a  Poisson  point 
process,  with  mean  rate  Event  2  also  occurs  according  to  a  Poisson  point 
process  with  rate  X2',  however,  some  of  these  latter  events  are  initiated  by  event 

1  while  others  occur  independently.  Suppose  that  event  2  is  initiated  by  event  1 
with  a  probability,  p,  and  at  a  random  time  delay  6.  Event  2  can  also  occur  at 
random,  with  mean  rate  of  occurrence  X2  as  before.  Then,  the  occurrence  of  event 

2  is  a  Poisson  process  (Wen,  1981),  with  mean  occurrence  rate  A,2  =  ^2  +  Pvi- 
Although  events  1  and  2  individually  are  Poisson  point  processes,  collectively,  the 
events  may  cluster  because  of  the  common  term  A,i. 

The  mean  rate  of  coincidences  of  such  clustered  events  (cf  Eqns  3.18)  is 
dependent  on  the  distributions  of  the  durations  of  each  event  and  on  the  delay 
time,©.  If  the  event  durations  and  6  are  modeled  by  exponential  distributions,  the 
mean  coincidence  rate  is,  approximately  (Wen,  1981), 

?112  =  ?'-2^l('C2  +  '^l)["i  +  Pvl/(^2^1  P-e]  (3.20) 

in  which  p©  is  the  mean  delay  time.  Note  that  the  first  part  of  Eqn  3.20  is  the 
same  as  Eqn  3.18c:  the  portion  in  brackets  represents  the  contribution  of 
stochastic  dependence  in  event  occurrences. 


29 


/\  -  initiia-ted  events 


Figure  3,2  Stochastically  dependent  event  occurrences 
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3.3  Illustration  of  Event  Combination  Analysis 

Statistical  characteristics  of  the  structural  loads  are  summarized  in  Table 
2  1  Type  II  loads  (Figure  2.1)  occur  intermittently  and  may  equal  zero  for 
significant  periods  of  Le.  If  the  pulse  process  is  intermittent,  the  probability 
that  the  process  is  nonzero  at  any  time  is  p=Xt.  For  a  combination  of  two 
intermittent  processes,  the  mean  rate  of  occurrence  and  the  duration  of  a 
coincidence  of  two  Type  II  loads  i  and  j  are  given  in  Eqns  3.18  and  3.19.  The  mean 
rates  of  occurrence  and  durations  of  loads  acting  individually  and  in  combination 
can  be  used  to  identify  those  combinations  of  events  that  lead  to  load 
combinations  requiring  codification  for  design.  The  principal  action-companion 
action  format  used  in  modern  limit  states  design  (cf  Eqns.  1.5  -  1.9)  requires  that 
each  of  the  time-varying  loads  assume,  in  turn,  the  position  of  the  principal  load. 

If  the  probability  of  a  joint  occurrence  of  a  companion  action  with  the  principa 
action  is  very  small,  then  the  companion  action  may  be  set  equal  to  zero  in  the 
load  combination  without  much  affecting  the  probability  that  the  design  load 
combination  will  be  exceeded.  Deterministic  load  combinations  involving  loads 
with  essentially  zero  probability  of  load  coincidence  do  not  enhance  structural 
reliability.  For  practical  purposes,  if  a  combined  load  event  has  a  coincidence 
probability  of  less  than  10-5/yr,  its  contribution  to  risk  is  negligible  and  there  is 
no  need  to  consider  it  in  a  load  combination  analysis. 


In  the  following  paragraphs,  the  occurrence  rates  of  several  possible 
combined  events  for  navigation  structures  are  illustrated.  Such  an  analysis  for  a 
particular  facility  may  enable  a  number  of  event  combinations  to  be  screened 
from  further  design  consideration. 

Oppratina  loarif^  end  wind.  The  mean  rate  of  occurrence  of  operating  loads  Hd 
and  Q  is  400/month  (4800/yr)  and  the  mean  duration  is  2  minutes  (3.8  x  10-6  yr). 
The  mean  rate  of  occurrence  of  significant  extratropical  winds  is  assumed  to  be 
2/yr  with  a  typical  duration  of  4  hr  (4.56  x  10-4  yr).  The  coincidence  of  operating 
load  and  wind  has  a  mean  rate  of  occurrence. 


;iow  =  (4800  X  2)  (3.8  x  10-6  +  4.56  x  10-4)  =  4.41 /yr 


TOW  =  (3.8  X  10-6)  (4.56  x  10-4)/(3.8  x  10-6  +  4.56  x  10-4)  =  3.77  x  10-6  yr 
pow  =  (4.41 /yr)  (3.77  x  10-6  yr)  =  1.66  x  10-5 


Onpratinn  loads  ^nd  earthQuak_e.  The  mean  rate  of  occurrence  of  potentially 
significant  earthquakes  (magnitude  5  Dr  greater)  is  assumed  to  be  0.02/yr  or  less 
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at  sites  in  the  upper  Mississippi,  with  a  mean  duration  of  30  seconds  (9.51  x  10-7 
yr).  The  coincidence  of  operating  load  and  earthquake  has  a  mean  occurrence  rate 

XoE  =  (4800  X  0.02)  (3.8  x  10-6  +  9.51  x  10-7)  =  4.56  x  10-4/yr 

Toe  =  7.61  x  10-7yr 

Poe  =  8.7  x  10-io 

Operating  loads  and  impact.  The  mean  rate  of  occurrence  of  significant 
impacts  at  miter  gates  is  about  0.19/yr,  with  a  duration  less  than  15  seconds 
(4.76  X  10-7  yr).  The  mean  occurrence  rate  of  the  combined  event  is, 

Xo\  =  (4800  X  0.19)  (3.8  x  10-6  +  4.76  x  10-7)  =  0.0039/yr 

Toi  =  4.23  X  10-7yr 

Poi  =  1.7  X  10-9 

This  analysis  assumes,  however,  that  an  impact  can  occur  when  the  gates  are  in 
motion;  if  lock  operation  precludes  such  an  event,  the  mean  rate  of  occurrence  of 
a  coincidence  would  be  zero. 

Earthquake  plus  impact.  Using  the  occurrence  statistics  from  previous 
examples,  the  mean  rate  of  occurrence  of  a  coincidence  of  these  events  is, 

Xe\  =  (0.02  X  0.19)  (9.51  X  10-7  +  4.76  x  10-7)  =  5.42  x  10-9/yr 

TEi  =  3.17  X  10-7yr 

Pei  =  4.3  x  10-5 

Earthquake  plus  flood.  Assume  that  the  mean  rate  of  occurrence  of 
significant  flooding  is  0.10/yr  and  with  a  mean  duration  of  2  days  ((0.0055  yr). 
The  mean  rate  of  occurrence  of  a  coincidence  of  extreme  flood  and  earthquake  is, 

;.ef  =  (0.02  X  0.10)  (9.51  X  10-7  +  0.0055)  =  1.10  x  10-5/yr 
Tef  =  9-51  X  10-7yr 

Pef  =  2.6  X  10-11 
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On  the  basis  of  the  above  illustrations,  only  the  combination  of  operating  load  and 
wind  would  appear  to  be  worthy  of  additional  consideration  in  developing  design 
load  combinations.  The  results  are  site-specific,  however,  and  additional  studies 
should  be  performed  before  drawing  any  definitive  conclusions. 
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4.  RECOMMENDATIONS  FOR  DESIGN 


While  the  event  combination  methodology  described  in  Section  3.2  might  be 
applied  in  specific  instances  for  evaluating  an  existing  structure,  for  general 
design  purposes  the  probabilistic  statement  must  be  transformed  to  a  set  of 
deterministic  criteria.  With  the  probability  distribution  of  Umax  determined,  a 
design  value,  Ud,  is  selected  with  an  acceptably  small  probability  of  being 
exceeded.  There  are  three  barriers  to  this  approach:  (1)  A  target  probability  for 
the  event  that  Umax  exceeds  Ud  must  be  established:  (2)  The  determination  of 
fractiles  of  a  distribution  is  difficult  and  incompatible  with  routine  design 
computations;  and  (3)  There  is  a  need  to  codify  the  design  process  to  apply  to  a 
wide  variety  of  design  situations  falling  within  the  scope  of  the  building  code. 

A  target  reliability  measure  (limit  state  probability  or  reliability  index) 
must  be  established  as  a  basis  for  design  and  evaluation  considering  the  events 
that  have  not  been  screened  out  by  the  process  in  Section  3.2. 

The  limit  state  probability  due  to  a  particular  event  for  a  component  or 
system  can  be  written  as 

P(F)  =  P(F|E)P(E)  (4.1) 

in  which  P(E)  =  probability  of  the  event  and  P(F|E)  =  limit  state  probability,  given 
the  event.  This  breakdown  of  the  limit  state  probability  into  its  constituent  parts 
is  useful  in  order  to  focus  attention  on  appropriate  strategies  for  damage 
prevention,  mitigation  and  control.  Reductions  in  failure  probability  and 
consequences  of  severe  events  can  be  accomplished  by  reducing  either  of  the  two 
probabilities  in  Eqn  4.1.  A  similar  approach  was  taken  in  developing  design 
criteria  to  reduce  the  risk  of  progressive  collapse  in  buildings  (Ellingwood  and 
Leyendecker,  1978)  or  to  reduce  risks  from  fires  (GIB  W1 4,1 983). 

The  event  probability,  P(E),  can  be  reduced  by  reducing  or  eliminating  the 
possibility  of  hazard  in  the  operation  of  the  facility.  This  can  be  done  by 
controlling  the  use  of  the  facility  or  by  training  and  education.  Generally,  it  is 
not  feasible  to  eliminate  all  sources  of  hazard.  Such  measures  to  reduce  P(E) 
generally  are  nonstructural  and,  indeed,  nontechnical  in  nature. 

The  probability  of  failure  given  a  structurally  significant  accident,  P(F|E) 
depends  primarily  on  the  nature  of  the  construction  of  the  facility.  Steps  taken 
by  the  structural  engineer  to  design  damage  tolerance  into  the  structure  through 
enhanced  strength  or  protective  systems  mainly  affect  the  probability,  P(F|E). 
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The  development  of  appropriate  structural  design  criteria,  therefore,  begins  with 
this  quantity. 

Event  probabilities  depend  on  the  nature  of  the  event.  For  normal  operating 
events,  P(E)  =  1.0  and  P(F)  =  P(F|E).  On  the  other  hand,  rare  or  accidental  events 
such  as  significant  barge  impact  occur  infrequently  (mean  rate  of  occurrence  on 
the  order  of  10"3  or  less,  cf  Table  2.1).  For  rare  or  accidental  events,  special 
measures  beyond  the  usual  ASD  or  LRFD  methods  (e.g.,  protective  or  sacrificial 
elements)  may  be  appropriate.  For  example  the  limit  state  probability  (on  an 
annual  basis)  for  significant  barge  impact  would  be, 

P(F)  =  10-3  P(FIE)  (^-2) 

To  maintain  economy  in  the  allocation  of  design  and  construction  resources, 
special  design  measures  should  be  selected  so  that  the  resulting  P(F|E)  is  such 
that  the  risk  of  failure  due  to  the  accidental  event  is  balanced  with  competing 
risks  from  other  operational  and  environmental  events. 

To  establish  a  frame  of  reference  for  these  risks,  the  limit  state 
probabilities  of  steel  or  reinforced  concrete  flexural  members  designed  to 
withstand  gravity  loads  are  on  the  order  of  0.0005  to  0.005  on  a  50-year  basis, 
(Ellingwood,  et  al,  1982;  Galambos,  et  al,  1982).  On  an  annual  basis,  then,  these 
limit  state  probabilities  are  of  the  order  10-5  to  10-^.  However,  the  flexural 
limit  state  usually  is  benign  and  not  life-threatening.  Because  of  the  severe 
economic  consequences  of  failure  of  a  navigation  facility,  it  is  reasonable  to 
require  P(F)  to  be  less  than  10-5.  |f  one  were  to  set  P(F|E)  =  0.01  in  the  above 
example,  the  annual  probability  of  structural  failure  would  be  on  the  order  of 
1 0-5,  placing  this  risk  in  the  low-probability  background  with  other  involuntary 
risks  (Wilson  and  Crouch,  1987). 

Criteria  for  risk-consistent  design  of  a  navigation  facility  should  be 
developed  to  be  consistent  with  a  set  of  performance  objectives  expressed  in 
probabilistic  terms.  A  conditional  limit  state  function,  g(  ),  must  be  postulated 
using  principles  of  mechanics  and  experimental  data  to  relate  the  load  and 
structural  response  variables  of  interest.  The  probability  distributions  of  each 
variable,  X,  must  be  determined.  The  (conditional)  limit  state  probability  is 
obtained  by  integrating  the  joint  probability  density  function  of  X  over  that  region 
where  g(  )  <  0.  In  the  example  above,  this  probability  would  be  compared  with  the 
target  limit  state  probability  P[F1E]  =  0.01;  if  the  computed  probability  is  greater 
than  the  target,  the  strength  of  the  structure  must  be  enhanced  in  some  way. 
Similarly,  design  requirements  corresponding  to  a  target  limit  state  probability 
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of  0.01  can  be  developed  (Galambos,  et  al,  1982;  Ellingwood,  et  al,  1982).  One 
can  consider  only  the  load  combination  aspect  of  the  safety  check,  leaving 
material  behavior  and  limit  state  considerations  to  responsible  material 
specification-writing  groups. 

Research  in  probabilistic  load  combination  analysis  has  shown  that  failures 
that  do  not  arise  from  negligence  usually  occur  when  one  of  the  time-varying 
loads  attains  its  maximum  value  while  the  other  loads  equal  their  point-in-time 
values  (Turkstra  and  Madsen,  1980).  This  suggests  that  Eqn  1.2  can  be  replaced 
for  practical  design  purposes  by  the  set  of  equations. 


Ud  =  Yd  D  -h  Yq  Qn  +  lYiQi  (4.3) 

The  load  combination  format  in  Eqn  4.3  is  referred  to  as  a  "companion  action" 
format.  The  factored  load  Yq  Q  is  denoted  the  principal  action,  while  the  terms  Yi 
Qi  are  the  companion  actions.  In  principle,  each  variable  load  must  be  positioned, 
in  turn,  in  the  position  of  the  principal  variable  load  in  order  to  determine  the 
maximum  combined  effect;  in  practice,  this  seldom  is  necessary.  Most  modern 
load  combination  schemes  invoke  a  companion  action  format  of  some  kind  (ASCE 
7-88,  1990;  Eurocode  No.  1,  1991).  The  development  of  specific  load 
combinations  is  outside  the  scope  of  this  report. 
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5.  CONCLUSIONS 


Event  combinations  for  design  or  evaluation  of  navigation  structures  should 
be  evaluated  probabilistically,  and  design  criteria  can  be  developed  to  be 
consistent  with  a  level  of  performance  expressed  in  probabilistic  terms. 

1.  The  probability  of  a  joint  occurrence  of  common  loads  at  their  peak 
values  is  very  small.  Thus,  when  considering  what  loads  should  be  combined,  a 
significant  reduction  in  the  combined  load  from  the  sum  of  nominal  design  values 
is  appropriate.  A  principal  action-companion  action  format  for  combining 
structural  actions  is  recommended  for  maintaining  uniform  reliability. 

2.  The  illustration  of  the  event  combination  methodology  in  section  3.3 
indicated  the  following: 

The  probability  of  a  coincidence  of  operating  loads  and  earthquake  is  on  the 

order  of  10-io. 

The  probability  of  a  coincidence  of  earthquake  and  impact  is  on  the  order  of 

10-15. 

The  probability  of  a  coincidence  of  earthquake  and  flood  is  on  the  order  of 

10-11. 


It  appears  as  if  combinations  of  these  events  are  very  rare  and  would  not  need  to 
be  considered  in  design  or  evaluation.  However,  further  study  of  these 
recommendations  would  be  desirable. 

3.  Load  combinations  for  dealing  with  accidental  events  should  be  based  on 
a  conditional  probability  of  about  0.01-0.1  that  the  combined  load  effect  is 
exceeded,  given  that  an  accident  occurs.  The  overall  unconditional  probability  of 
structural  failure  using  this  design  approach  is  estimated  to  be  on  the  order  of  10- 
5/yr  or  less. 

4.  Research  should  be  undertaken  to  develop  a  protocol  for  combining 
extreme  events  in  design  of  navigation  facilities. 
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